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Abstract 

This paper presents a new distributed control framework to coordinate inverter-interfaced distributed energy resources (DERs) 
in island microgrids. We show that under bounded load uncertainties, the proposed control method can steer the microgrid to 
a desired steady state with synchronized inverter frequency across the network and proportional sharing of both active and 
reactive powers among the inverters. We also show that such convergence can be achieved while respecting constraints on 
voltage magnitude and branch angle differences. The controller is robust under various contingency scenarios, including loss of 
communication links and failures of DERs. The proposed controller is applicable to lossy mesh microgrids with heterogeneous 
R/X distribution lines and reasonable parameter variations. Simulations based on various microgrid operation scenarios are 
also provided to show the effectiveness of the proposed control method. 
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1 Introduction 

Microgrids are low voltage power networks comprised 
of distributed generations (DGs), energy storages sys¬ 
tems (ESSs), and loads that can operate in either grid- 
connected or island mode. Distributed generation con¬ 
tributes on-site and clean energy, which is expected to 
make power networks more robust, efficient and envi¬ 
ronmentally friendly [1,2]. Energy storage systems are 
considered as an important resource to benefit the power 
networks by smoothing real time imbalance between gen¬ 
eration and demand [3]. Some storage devices such as 
freewheel and battery packs can be integrated with in¬ 
termittent DGs to regulate the power injection to a 
power network [4, 5]. Demand side appliances such as 
plug-in hybrid electric vehicles (PHEV) and thermostat¬ 
ically controlled loads (TCLs) can also be viewed as en¬ 
ergy storage resources. Those “storage” appliances can 
be coordinated to provide ancillary services to the main 
grid [6-8]. The proximity of DGs and ESSs to loads in 
a microgrid allows for a transition to the island mode 
during faults on the main grid. Such a transition may 
also be triggered by efficiency or reliability incentives, 
(see [9,10]). 
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Distributed energy resources (DERs) such as DGs and 
ESSs connect to the microgrid through DC/AC or 
AC/AC inverters. During the island mode, the invert¬ 
ers are typically operated as voltage source inverters 
(VSIs). These VSIs need to be controlled cooperatively 
to achieve desired performance and reliability proper¬ 
ties. In AC networks, voltage magnitude and angle dif¬ 
ference between connected buses should be regulated in 
some bounded ranges for system security and stability. 
Frequency synchronization to a nominal value is also 
crucial for grid connection and stability purposes. Be¬ 
sides frequency and voltage regulation, sharing of active 
and reactive power is also considered as important con¬ 
trol objectives in microgrids [11,12]. They require that 
the power injection into the microgrid from DERs is 
proportional to the nominal value defined by economics 
or other incentives, while satisfying load demands [13]. 
Power sharing enables effective utilization of limited 
generation resources and prevents overloading [14]. 

To achieve the aforementioned objectives, a microgrid is 
typically controlled using a hierarchal structure includ¬ 
ing primary, secondary, and tertiary controls [15-18], 
which is similar to the one used in the traditional power 
systems. The primary droop control of a microgrid main¬ 
tains the voltage and frequency stability while balanc¬ 
ing the generation and load with proper power sharing. 



The secondary controller compensates the voltage and 
frequency deviations from their reference values. The 
tertiary control establishes the optimal power sharing 
between inverters in both islanding and grid-connected 
modes. 

The primary droop is generally a decentralized controller 
that adjusts the voltage frequency and magnitude of 
each inverter in response to active and reactive power de¬ 
viations from their nominal values. Various droop meth¬ 
ods are proposed to achieve proportional active and re¬ 
active load power sharing [11,19-24]. However, this is 
often achieved at the cost of sacrificing other control ob¬ 
jectives such as voltage and frequency regulation. The 
secondary control utilizes either centralized or decentral¬ 
ized communication infrastructures to restore frequency 
and voltage deviation induced by the primary droop. 
Most of the existing secondary control methods require 
centralized communications [25-27]. On the other hand, 
decentralized secondary control has recently been pro¬ 
posed to avoid single point of failure [28]. The combined 
operations of the primary and secondary control require 
separation of time scale, resulting in slow dynamics that 
cannot effectively handle fast-varying loads [29]. In ad¬ 
dition, the secondary control may destroy the propor¬ 
tional power sharing established in the primary control 
layer [30]. One possible solution is to adopt distributed 
or decentralized control structure for primary and sec¬ 
ondary control layers to improve performance and sup¬ 
port plug-and-play operation of the microgrid [18]. 

Many existing primary and secondary control methods 
rely on small signal linearization for stability analysis, 
which is vulnerable to parameter variations and change 
of operating points. Only several recent works [13,30,31] 
have rigorously analyzed the stability of microgrid with 
droop-controlled inverters. In particular, [30] derives a 
necessary and sufficient condition for the stability under 
primary droop control. The authors have also proposed a 
distributed averaging controller to fix the time scale sep¬ 
aration issue between the primary and secondary control 
layers. In [13] and [31], stability conditions of lossless 
mesh microgrids have been provided. Despite their ad¬ 
vantages, these nonlinear methods still suffer from sev¬ 
eral common limitations. First, all the nonlinear anal¬ 
yses mentioned above only focus on lossless microgrids 
with purely inductive distribution lines. The results may 
not be applicable for microgrids with heterogeneous and 
mixed R/X ratio lines, which is common in low voltage 
microgrids [25]. Secondly, since only frequency droop is 
carefully analyzed, reactive power sharing is often not 
guaranteed. 

To address the aforementioned limitations of the exist¬ 
ing works, we propose a distributed control framework 
to coordinate VSIs in an island AC microgrid. The pro¬ 
posed control adjusts each inverter frequency and volt¬ 
age magnitude based on the active/reactive power mea¬ 
surements of its neighbors. We first show that the par¬ 


ticular control structure ensures that any equilibrium 
of the closed-loop system results in the desired power 
sharing and frequency synchronization. Secondly, condi¬ 
tions for power sharing and synchronized frequency re¬ 
specting voltage constraints are provided. The proposed 
controller can be applied to both radial and mesh mi¬ 
crogrids with mixed R/X ratios. Furthermore, the pro¬ 
posed controller requires no separation of time scale and 
can tolerate reasonable parameter variations. To the au¬ 
thors’ knowledge, most existing control framework can¬ 
not achieves active/reactive power sharing while respect¬ 
ing voltage and frequency regulation for a mesh micogrid 
with mixed R/X ratio lines. 

To demonstrate the robustness of the proposed dis¬ 
tributed controller, we also study the control perfor¬ 
mance under partial communication failures and the 
plug-and-play operations. We will show that as long 
as the communication network remains connected, all 
the desired properties including power sharing and 
frequency and voltage regulation still hold in these con¬ 
tingency scenarios. This effectively demonstrates the 
robustness of the proposed distributed controller. It is 
worth to mention that the proposed framework may re¬ 
quire faster communications among the VSIs than the 
traditional secondary control. However, such commu¬ 
nication requirement is reasonable for most microgrid 
control systems [32-34]. 

The rest of this paper is organized as follows. Section 2 
formulates the microgrid control problem. Sufficient con¬ 
ditions for the solvability of the proportional power shar¬ 
ing problem respecting voltage constraints are also pro¬ 
vided. The proposed distributed control framework is 
developed in Section 3. Robustness of the distributed 
controller under loss of communication links or failures 
of DERs is studied in Section 4. In Section 5, we vali¬ 
date the proposed controller through simulations under 
various microgrid operating scenarios, including abrupt 
changes of loads and loss of one VSI. Some concluding 
remarks are given in Section 6. 

Notation Define R + and R_ as positive and negative 
real numbers, respectively. Denote [n] := {1,2,..., n}. 
Given a set V, let |V| and 2 V be its cardinality and power 
set, respectively. Denote the diagonal matrix of a vector 
x as diag(x). For a set of vectors Xi, i G X, let {xi, i € 1} 
be the augmented vector of Xi collecting all i G X. Given 
a polyhedron B G R”, let v(B) be the vertex set of B. For 
a closed set F C R", int/F) and dF are the interior and 
the boundary of F. The distance between a point_/ G R™ 
and the set F is denoted as d(f, F) := inf{ | ]/ — /] | 2 ]/ G 
F}. Define l n G R" and 0 n G R" as the vectors with 
all the elements being ones and zeros, respectively. For a 
symmetric matrix A , let A(A) and A(A) be the spectrum 
and minimal eigenvalue of A, respectively. Denote A®B 
as the tensor product between matrices A and B. Let 
null(A) be the null space a matrix A. 
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Fig. 1. A connected microgrid network. The dash lines rep¬ 
resent the communication links and the solid lines represent 
the distribution lines connecting VSIs, loads, and the main 
grid. 


Xi can be fully controlled. In contrast to Xi, voltage at 
load buses Xl is uncertain. The voltage dynamics of the 
load buses are assumed to satisfy the following condition 

Assumption 1 | |cc^1 1 2 < /c||a;j|| 2 , for some n > 0 

where k > 0 is a constant determined by load prop¬ 
erties and microgrid topology. More explanation about 
this assumption is provided in Appendix A. Throughout 
this work, we will focus on constant power or constant 
impedance loads so that Assumption 1 holds. 

Remark 1 Under multi-agent or centralized control 
framework, load voltage is often assumed to be measur¬ 
able and known during controller design [39,40]. In this 
paper, we consider a more general scenario, where load 
voltage is viewed as an unknown variables with only 
Assumption 1 being involved in the controller design. 


2 Problem Formulation 

In this paper, we consider a connected island microgrid 
network as shown in Fig. 1. An island microgrid is repre¬ 
sented by a connected and undirected graph Q = (V, £), 
where V is the set of buses (nodes) and £ C V x V 
is the set of distribution lines (edges) connecting the 
buses. The set of buses is partitioned into two parts, 
inverter buses V/ and load buses Vj> Let raj = |V/|, 
riL = \Vl\ and n = |V|. The magnitude and phase an¬ 
gle of the bus voltage are denoted as E t and 0i, respec¬ 
tively. Let Xi = [6i,Ei\ T be the state vector at bus i, 
and let xi = {xi,i £ V/} and xl = {xi,i £ Vj,} be 
the inverter bus state vector and load bus state vector, 
respectively. The overall system state vector is denoted 
by x = [xJ,x^] T and will be referred to as the system 
voltage profile. For each bus i £ V, let Pi and Qi be the 
active and reactive power injections at bus i. Given the 
admittance matrix Y £ C nx ” of the microgrid, the 
active and reactive power injections are related to the 
voltage profile x by the power flow equations [35] 

f Pi{x ) = Ei YijEjCos(9i — 0j — (j>ij) 

0 i ( X ) — P'i y ' EjSlll{fii 0 j ), 


where Yjy = ||Yyj| 2 and <pij = ZY .y are the magnitude 
and the phase angle of the admittance matrix element 


We distinguish the voltage at inverter and load buses in 
our formulation due to their different characteristics. For 
inverter buses, there are standard methods to control 
the voltage magnitude and frequency ( [36], [37]). Typ¬ 
ically, these methods can track a given inverter voltage 
reference almost instantaneously. Therefore, an inverter 
is often modeled as a controlled voltage source behind a 
reactance [38]. We also adopt such a model and consider 


Given nominal active and reactive power injections P* 
and Q*, i £ Vj, it is desired that the power injection 
of the inverters share uncertain loads proportionally to 
their nominal value: 


Definition 1 (Proportional Power Sharing) The 

active and reactive power are proportionally shared 
among the buses j £ V/ if 



Pk{x) Qj(x) 
p; ’ Qj 


Qk (•£) 

~qT 


j,k £ Vi. 


( 2 ) 


The power sharing condition (2) imposes a constraint on 
the system voltage profile. We define this constraint set 
as 


x £ Xs := {x|Eq. (2) holds}. (3) 

In addition to condition (3), the control and operation 
of a microgrid has to respect its branch angle difference 
limits and voltage magnitude constraint. The branch an¬ 
gle difference between all connected buses is typically 
required to be bounded by a given constant 7 £ [ 0 , f). 
The upper bound 7 is derived based on the maximum 
current allowable on each distribution line I^ ax (see [41] 
and [42]). In addition to branch angle difference, the 
voltage magnitude also needs to stay inside some secure 
operation range [35]. Denote [E^Ei] as a given desired 
range of the voltage magnitude of bus i. Both branch an¬ 
gle difference and voltage magnitude requirements im¬ 
pose a constraint on the system voltage profile defined 
below: 


X 0 = {x\\d z - Ofi < 7 , V{f,j} £f}, 
X E = {x\E_i <E l <E i Vi£ V}, 


These two constraints will be referred to as the security 
constraint of the microgrid 
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Definition 2 (Security Constraints) We say that a 
microgrid satisfies the security constraints if 

xeX c -.= x E nx e . (4) 


In addition to proportional power sharing and voltage 
regulation, another important nricrogrid control objec¬ 
tive is known as frequency regulation. Frequency regu¬ 
lation is defined as synchronization without deviations 
from the nominal value, specifically, 9i = u>o, Vi € V, 
where coq is a predefined nominal frequency of the mi¬ 
crogrid. Note that the power flow equations (1) and the 
security constraints are invariant with respect to rigid 
rotation of 9i of all buses. We can select a reference 
frame rotating at angular frequency u >o while preserving 
all properties in Eq. (l)-(4). With the rotating reference, 
the frequency regulation condition is reduced to 

9i = 0, Vi G V. (5) 


As discussed in [17], the controller design to meet the 
requirements specified in Eq. (3)-(5) inevitably requires 
communication networks. In this paper, we employ a 
distributed communication structure similar to [ 12 ] 
and [29], where each inverter can communicate with its 
neighboring inverters to share its local measurements as 
shown in Fig. 1. Let Q c = (V/, £ c ) be a connected simple 
graph of the communication network, where each edge 
{i,j} G £ c represents an available communication link 
between buses i and j. Let A /) := {j\{i,j} G £ c j U {?'} 
be the set of neighbors of bus i, (including bus i itself). 
An inverter i has access to the measurements at every 
inverter bus j G A), including Pj and Qj . 

Since each inverter is modeled as a VSI, the microgrid- 
level coordination control for each inverter i reduces 
to the determination of appropriate voltage frequency 
and magnitude setpoints. The actual frequency 6, and 
magnitude Ei can track these setpoints almost instan¬ 
taneously. The challenge here lies in that the load is 
uncertain and different load conditions require different 
voltage profile xj in order to satisfy constraints (3)-(5). 
Define S, = {Pt/P* ,Qi/Q*] T , = {Sj,j G A /]} and 

E jy/v = {Ej,j G Ahi}. Our goal is thus to design a con¬ 
troller for each inverter i that can automatically find 
the desired voltage vector x, based on local information 
(Sx, • E j \r i ). Towards this end, we propose to dynami¬ 
cally update Xi as follows 

±i = m {Sm, (a:), E ^), Vi G V/ 

( 6 ) 

subj. to x(t) G X c , Vt > 0, 

where Hi is the control law of inverter i to be designed. 
Note that the above control structure corresponds to 


directly assigning frequency 9i based on local informa¬ 
tion, while dynamically updating voltage magnitude Ei 
through simple integrator dynamics. Such structure is 
commonly used in the literature of microgrid control, 
(see [12] and [43]). The constraint x(t) G X c is imposed 
to ensure that the security constraints are always satis¬ 
fied. 

Define X e := X c n X$, then conditions (3) and (4) hold 
when x G X e . Under Assumption 1, condition (5) holds 
when ±i = 0 for all i G V/. Our goal becomes to design¬ 
ing Hi such that X e forms an equilibrium set of system 
(6). In addition, we also want to achieve an exponential 
convergence to X e for some initial a:(0). If such a con¬ 
troller is found, it can steer the microgrid to the desired 
steady state where conditions (3)-(5) hold. In the rest of 
this paper, we will first develop the control law Hi such 
that X e forms the equilibrium set of the system ( 6 ), and 
then derive conditions to ensure the exponential conver¬ 
gence of X e . 

3 A Distributed Microgrid Control Framework 

In this section, we propose a distributed control frame¬ 
work to coordinate the inverters in an island AC mi¬ 
crogrid to accomplish the control objectives (3)-(5). We 
first provide sufficient conditions to ensure X e 0. A 
control design framework is then developed. 

3.1 Existence of Solutions 

A minimum requirement for the controller design is the 
existence of a voltage profile x satisfying all the con¬ 
straints, i.e., X e 0. Existing methods in the literature 
often directly assume this condition holds ( [ 20 ], [ 22 ]). 
Here, we provide a brief discussion and a set of sufficient 
conditions to guarantee the non-emptiness of X e . The 
existence of the voltage profile x satisfying conditions 
(3) and (4) involves solving nonlinear algebraic power 
flow equations (1). We revisit a classical result in the 
following. 

Lemma 1 [41] Suppose that the following conditions 

hold 

(a) The microgrid is connected, 

(b) The admittance matrix Y is symmetric, 

(c) 2 Ej > Ek for all j , k G V, 

(d) /jr x < f B jk for all {j , k} G £, 

( e ) Ej ( — Bjj + Sfcev 7 Bjk) > ( J2kev\{j} BjkEk .) for 
all j G Vl and the strict inequality holds for at least 
one j G Vl, 

(f) Pi G [Ei E] Vi G V, Qi G [Q. Qi \ Vi G V L , 

where Bjk = Yjksm(<j>jk) and E%,Ei , Q.,Qi are con¬ 
stants determined by microgrid parameters including line 
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impedance and bounds of voltage regulation. Then there 
exists a solution to Eg. (1) such that x € X c . 

Lemma 1 is in fact a direct consequence of Theorem 4 
in [41]. Readers are referred to [41] for the proof and 
details of finding P 2 , Pi, Q . and Q,;. 

Remark 2 For simplicity, several conditions in [41] re¬ 
lated to loads serviceability are not included in Lemma 1. 
Since a transition to island mode is enabled only when 
the DERs can provide sufficient power to loads in the 
microgrid, the serviceability requirement is satisfied in¬ 
trinsically for this work. 

Remark 3 If the nominal active power injection at in¬ 
verter buses satisfying P* G [P t Pi], \/i G V/, Lemma 1 
implies X e = X$ H X c ^ 0. Since we focus on microgrid 
control problem with given P* and Q* , we thus assume 
that P* and Q* selected in the tertiary control layer are 
chosen such that P* G [P, ; P,], \/i G Vj. In this way, 
X e is nonempty if conditions in Lemma 1 hold. We can 
then focus on designing controller to steer the microgrid 
to X e . 

3.2 Distributed Controller Design 


Define O = span{u p ,u g }, where v p ,v q G R 2ni , v p = 
[1, 0,1, • • • , 0] T and v q = [0,1, 0, • • • , 1] T . The following 
proposition shows that under some mild conditions, ev¬ 
ery equilibrium point of system ( 8 ) satisfies the control 
objectives (3)-(5). 

Proposition 1 If Assumption 1 holds and null(K )C O, 
the following statements are equivalent 

(a) The microgrid with dynamics (8) is in steady state 
where Xi =0. 

(b) The desired conditions (3)-(5) hold. 


Proof. Since null(L) = l ni , null(Z)= O. Null(/vL) = 
O follows directly from null(/V)C O and null(L)= O 
With null( ATL)= O , we have 

±i = 0 such that x G X c (9) 

<t=> x G {a;|,S7(;r) G O} n X c 4$ x G X e , 

The equivalence between statements (a) and (b) follows 
from Eq. (9). □ 


We start our controller design from a simple property of 
a connected graph Q c . Let L G be the Laplacian 

of Q c . The null space of L is span{l nj } because Q c is 
connected. Observing that l n/ has a close relation with 
condition (3), we design /q(-) as a simple linear feedback 
in terms of <SV; i n the following form 

f ii(t) = KiJ2 jeVl L(i,j)Sj(x(t)) 

\ subj. to x{t) G X c , Vf > 0, 

where Ki G R 2x2 is the local control gain matrix at 
bus i to be designed. Define Sj = {Si, i G V/}, and 
S = {S',;, i G V}. Let K = diag{I\i,i G V/}, and L = 
L®l 2 - The dynamical model of the microgrid under the 
proposed inverter control (7) becomes 


ii(t) = KLSi(x(t)) 
subj. to x(t) G X c , Vf > 0, 


Remark 4 The proposed control structure ( 8 ) does not 
depend on the voltage magnitude information Ejy. that 
is also available at bus i. We will show later that such 
a control structure is already sufficient to ensure con¬ 
vergence to X e . In principle, the magnitude information 
Ejy t can be used to further improve the control perfor¬ 
mance, especially for voltage regulation. However, we 
will not study such extension in this paper. 


Proposition 1 reduces the microgrid control problem 
with numerous requirements to the study of exponential 
convergence to X e of system ( 8 ). We will therefore focus 
on analyzing system ( 8 ). 

3.3 Analysis of System (8) 

In this subsection, we derive the conditions of exponen¬ 
tial convergence to X e where all the desired conditions 
(3)-(5) follow. The exponential convergence to T e of sys¬ 
tem ( 8 ) is challenging in general due to the nonlinearity 
of the underlying system and the uncertainty of the load 
bus states Xl ■ Instead of directly analyzing system ( 8 ), 
we apply the chain rule to obtain the dynamics of Sj 
under the proposed control strategy 

f S^xft)) = J^^KLS^x^)) + J L , x ( t )X L {t) 

\ x{t) G X c , Vf > 0, 

where Jj^ x and J^ x are the Jacobian matrices of Si(-) 
evaluated at x with respect to Xi and Xl , respectively. 
Notice that at every time instant and for any xl, Eq. (10) 
describes the dynamics of Sj when the dynamics of Xj is 
given by Eq. ( 8 ). According to Proposition 1, the conver¬ 
gence of Si(x(t)) to O of system (10) implies the conver¬ 
gence of the state trajectory x{t) to X e . The close rela¬ 
tion between these two stability properties motivates us 
to focus on system (10). To simplify notation, we define 
z(t) = S T {x(t)), B(t) = J I>x ( t ), and J = {Ji, x \x G X c }. 
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System (10) can then be written as a linear time varying 
(LTV) system 

z(t) = B(t)KLz(t) + w(t ) 

B(t) £ J, 

where w(t) = JL,x(t)^L(t) is considered as a disturbance 
of system (11). With this notation, system (11) becomes 
a stand alone dynamic system with state variable z sub¬ 
ject to unknown disturbance w(t). Note that in system 
(11), w is quadratically bounded by d(z, O) 

IMh = \\Jl,xXl\\2 < K||Ji,x||2||a:/||2 
= k\\JlMKLz\\ 2 < C\\Lz\\ 2 
= C -d(z,0), 

where £ G K+ is a constant depending on system param¬ 
eters as well as control gain K. Robust stability of the 
equilibriums of systems with bounded noise was studied 
in [44], which is reviewed in the following 

Definition 3 The set O is robustly stable of system (11) 
with degree £ if O is globally exponentially stable for all 
w such that ||w ||2 < C d(z, O). 

To analyze robust stability for O of system (11), we apply 
a standard change of coordinates. Define a change of 
coordinate matrix equation 

T= [vi,v 2 ,..,V2n-2,-n\-, M ], (12) 

IM 2 IKII 2 

where the first 2n — 2 vectors are arbitrary vectors such 
that T is an orthogonal matrix. Let z = T~ x z be the 
state vector in the new coordinate system. The LTV 
system (11) becomes 


reduced order system of (13) 

z(t) = A 11 (t)z(t)+w(t) 

-dll (t) £ A, 

where A := {G(B)\ B £ J}. Similar to w in system 
( 11 ), w is quadratically bounded by the state z shown 
in the following 

IMI2 < IHI2 < Cll^lb- 

The following lemma shows that exponential conver¬ 
gence to X e of system ( 8 ) follows if the origin of system 
(15) is robustly stable. 

Lemma 2 If the origin of system (15) is robustly stable 
with degree <(, then there exists an non-empty X cs C X c 
such that for all x(Q) G X CiS ,a:(t) exponentially converges 
to X e for system (8). 


Proof. Since systems (15) and (13) share the same dy¬ 
namics in the space K 2n/ \ O , robust stability of the ori¬ 
gin of system (15) implies O is robustly stable of sys¬ 
tem (13) with degree f. In addition, robust stability of 
O of system (11) (or system (13) ) guarantees the tra¬ 
jectory of z(t) is bounded. Define Sj = {<S7(:r)|a; £ X c }. 
With the bounded trajectory of z(t) = Si(x(t)), there 
exists X CtS C X c such that for all x(0) £ X CtS , Sj(x(t)) G 
S/ Vf, which implies x(t) £ X c \/t. Therefore, for all ini¬ 
tial x(0) G X c>s , z(t) = Sj(x(t )) converges to O with 
x(t) £ X c for all time in system (11) if system (15) is ro¬ 
bustly stable. Recall that Eq. (11) describes the dynam¬ 
ics of Sj(-) when the dynamics of Xj is given by Eq. ( 8 ). 
We then conclude that for all x(0) G V c>s , x(-) expo¬ 
nentially converges to X e for system ( 8 ) if the origin of 
system (15) is robustly stable. □ 


z(t) = T~ 1 B(t)KLTz(t) + T~ 1 w(t) 
B(t) £ J. 


We now provide a set of sufficient conditions for robust 
stability of the origin of system (15) 


Since the last two coordinates of the new basis V span 
O , the last two column vectors of LT are zeros and 


Proposition 2 The origin of system (15) is robustly sta¬ 
ble with degree 1 ) if there exist £ R+, U = U T >- 0 
such that 


T~ 1 B(t)KLT = 

An(t) 0 

, (14) 

'AhU + uAn +e(I + fU u " 


A 2 i(t) o_ 


U -el 


-< 0 , 


(16) 


where An (•) G R2(n,-l)x2(n7-l) and i 2) (.) g R 2x2(n J -l) 
Considering that the dynamics of d(z(-),G) of system 
( 11 ) is irrelevant to the last two coordinates of the 
state z of system (13), we focus on a reduced order 
system of (13) with state vector z = [I 0]z £ 

Dehne w = [I 0]T _ 1 r/; G IR 2 (™ 7-1 ) and G(B ) = 
[I 0 ]T~ 1 BKLT[I 0] T G K2(n / -l)x2(n J -l) ) we have a 


for all An G A. 


Proof. The proof is similar to linear time invariant sys¬ 
tem case discussed in [44]. Eq. (16) can be derived by 
quadratic Lyapunov function argument. If there exist a 
Lyapunov function V(z) = z T Uz such that for all An £ 
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A, V{z) < —£V(z), then the origin of system (15) is ex¬ 
ponentially stable. The conditions for V{z) < — fV(z), 
Viu G A such that ||io ||2 < Cl 1^1 b are shown in the 
following 


V(z) < -£V s.t. ||w|| 2 < CPIh 


" -T " 

Z 


A^U + U An 

+ fU u 


Z 

- T 
W 


U 

0_ 


W 



-a o' 


z 


0 1 


U) 


" ~T" 

Z 

\A^U + UAu+etl + tU U 1 

Z 

- T 

W 

1- 

c: 

1 

CT\ 

1 ^ 

W 


< 0 . 


for all An G A. S-procedure is applied for the last step, 
which completes the proof. □ 


Note that if An is polytopic, the condition (16) can be 
formulated into bilinear matrix inequalities (BMIs). The 
condition can then be checked numerically. However, An 
is not polytopic, so we will instead develop a way to find 
a convex set containing An in the next subsection. 

With Lemma 2 and Proposition 2, conditions of expo¬ 
nential convergence to A e of (8) can be obtained: 

Theorem 1 If Assumption 1, hypotheses in Lemma 1 
and Eq. (16) hold, then there exists an non-empty X cs C 
X c such that for x(0) G X CtS , the microgrid (8) converges 
exponentially to the set X e where the control objectives 
including (3)-(5) are all satisfied. 


Proof. The origin of system (15) is robustly stable due 
to Proposition 2. By Lemma 2, robust stability of system 
(15) implies the existence of X c ^ s C X c such that for all 
a:(0) G A CjS , the trajectories converge to X e of system 
(8). Since Eq. (16) ensures null(A")C O, X e is equivalent 
to desired control objectives (3)-(5) from Proposition 1 

□ 


3-4 Feedback Gain Design 

In this subsection, we propose a constructive way to find 
a feedback gain K satisfying Eq. (16). The difficulty lies 
in checking the feasibility Eq. (16). As discussed in the 
last subsection, the robust stability condition in Eq. (16) 
can not be directly formulated into BMIs because A is 
not polytopic. We will first derive a convex hull con¬ 
taining A by analyzing the Jacobian of the power flow 
equations (1) so that Eq. (16) can be checked by solv¬ 
ing several BMIs. Secondly, instead of only checking the 
feasibility, we formulate the BMIs into an optimization 
problem to enhance the robustness. 

Define C, C £ and Xc i as a polyhedron replacing the 
constraint 1 0i — 9j\ <7 in X c by 6i — 9j =0 for all {i,j} G 
Ci. Let $ := | + 7 < §}. An approximated 

convex hull can be found by the following Proposition. 

Proposition 3 If every entry of the admittance matrix 
Y satisfies faj G {±f, $}, the upper and lower bounds 
of every entry (i,j) of Ji, x G 3 are 

) = ma x{Ji tZ (i,j),z G Z}, (17) 

i Li tX (i,j) = min {J ItZ (i,j),z G Z}, 

where Z = \J Ci&2 £ v ( x cd- 


Proof. According to the power flow equation (1), the 
derivative of the active power injection at bus i with 
respect to different variants are 


dPfx) 

d9i 

dPi(x) 

dEi 

dPi(x) 

dOj 

dPi(x) 
dEj 


Ei ^ ) Vj y'Ay sin(YA, 0j fij ) • 

J 6 V\* 

2 A, Vj, cos( 0,r ) + ^ ^ YijEjCos(6i 0j fiij) ■ 
jev\i 

EiYijEj sin( 0, 0j (),j ), j f 1 i, 

EiYij cos(6»i 0j f ij ), j i. 

(18) 


The result of Theorem 1 is robust with respect to small 
variations of system parameters. As long as the pertur¬ 
bations of the admittance matrix are small enough such 
that Ji tX G 3, the exponential convergence to X e for 
some x(0) G X CiS still follows from Theorem 1. Different 
from most of the literature, the proposed controller can 
be applied to mixed R/X ratio distribution lines and 
general microgrid topology including acyclic and mesh 
networks. Furthermore, the controller can meet all the 
main control objectives without the separation of time 
scale, which distinguishes it from the mainstream droop 
control methods. 


Since for all entries in the admittance matrix Y, <pij G 
{±((, < I ) }, every summand of the Jacobian elements in 
(18) has the maximal and minimal points at x G Z if x G 
X c . In addition, every summand corresponds to different 
set of Ej and 9i — 9j in each Jacobian element in (18), 
the Jacobian elements have the maximal and minimal 
points at x G Z if x G X c . The same conclusion for the 
reactive part is reached by a similar argument. □ 


Given the upper and lower bounds of every entry of 
Jj x G 3, a convex hull containing 3 can be found. 
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Let F be a set {—1,1}. Define Qi £ p' 2niX2ni such that 
= F 2nix2 " J , where l = 2 4 < Define D Si £ 
R :2n,x2ni such that 

D ei {j, k ) = Ji }X (j, k ) if Qi(j, k) = 1 , 

D ei (j , k) = J IiX {j, k ) if Qi(j, k) = - 1 . 

The following lemma is a simple consequence of Propo¬ 
sition 3 and provides a convex hull containing J. 

Lemma 3 If every entry of the admittance matrix Y 
satisfy <fij € {±§, $}, the convex hull 

V:=CO{D ei ,...,D g J 

contains J. 

The results of Lemma 3 allows us to replace Eq. (16) 
in Theorem 1 by BMIs. Instead of only finding K such 
that Eq. (16) is feasible, we propose to design K by 
solving the following optimization problem subject to 
BMI constraints 

maximize £ 

K > u >*£ (19) 

subject to Mi(e, £,£) V 0 , 

where e, £ € R + and 

in,i = [I,0]T~ 1 D ei KLT[I,0] T , 

^ _ Ai lti U + UA llt i + efl + fU U 

U -el\ ’ 

U = U T y 0, i £ [/]. 


Remark 5 The maximization of £ is for the purpose of 
improving the convergence rate. Notice that maximiz¬ 
ing £ is only meaningful if an upper bound of 11 YY 112 is 
imposed. More discussions on the upper bound will be 
included in the next section. 

The following corollary is a direct consequence of 
Eq. (19) and Theorem 1 

Corollary 1 If Assumption 1, the hypotheses in Lemma 
1 hold, and the optimal solution K for Eq. (19) is im¬ 
plemented in Eq. (8), then there exists an non-empty 
X CtS Q Y c such that for x(0) € X C:S , the microgrid (8) 
converges exponentially to the set X e where the control 
objectives including (3)-(5) are all satisfied. 

We want to comment that the convex optimization prob¬ 
lem (19) subject to BMI constraints is NP-liard to solve 
in general. However, efficient algorithms [45] and [46] are 
available if an initial feasible solution can be found. We 


can first substitute U by some simple positive definite 
matrices as an initial guess and check the feasibility. If U 
is a feasible solution of Eq. (19), algorithms [45] or [46] 
can be applied to find a local optimal solution. Those 
algorithms only involve several linear matrix inequali¬ 
ties (LMIs) instead of BMIs, which can be effectively 
solved by using various existing convex optimization al¬ 
gorithms [47,48]. 

Remark 6 (Trade-off between Complexity and Robust■- 
ness) Let J := CO{Ji. z ,z £ Z}. The convex hull J 
does not necessarily contain J. However, there exists a 
compact set X c ^ s C X c such that J D {Ji, x ,x £ X C}S }. 
Hence, if one solves Eq. (19) by replacing v(T>) with 
Corollary 1 still implies exponential convergence 
to a subset_of X e for some initial x(0). The benefit of re¬ 
placing v(T>) with v(J) is reduction of complexity, where 
\v(f)) \ ~ 2 4 ”? while |u(jf)| ~ 2"-3^. The difference be¬ 
comes evident for large n/ and one may prefer to solve 
the_optimization problem (19) by substituting v(T>) to 
v(S) for microgrids with larger nj. 

As discussed in Remark 6 , the number of the BMI con¬ 
straints in the optimization problem (19) grows expo¬ 
nentially with respect to ?ij. Finding K for microgrids 
with large number of bus can become challenging even 
one replacing v(T>) by v{J). Here, we observe that ma¬ 
trices in J are block diagonal. The number of blocks 
equals to the number of groups of inverter bus separated 
by load bus. For example, all matrices of J of IEEE 14 
bus system shown in Fig. 2 has three diagonal blocks. 
One of the blocks corresponds to 1 th — 3 rd inverter buses, 
and the values of the entries depend on the states of 
1 th — 5 th buses, which are covered with the same color in 
Fig. 2. Similar argument applies to the other two blocks. 

Following the motivating example above, we introduce 
several notations to define properties of every block in 
J. Let Ci be the set of inverter buses associated with 
block i. The dimension of each block is then given by | c* |. 
Similar arguments to Proposition 3 and Lemma 3 are 
made to define the convex hull Vi associated with block 
i. The number of vertices of 2?,; is then f = 2 4 I C; I . Notice 
that the large number of BMIs in Eq. (19) is originated 
from the number of vertices of the convex hull V. In the 
following lemma, we will show that each block in the 
block diagonal Jacobian matrix can be viewed separately 
in finding a feasible solution of Eq. (19), resulting in 
much less number of BMIs involved. 

Lemma 4 Every I\ that satisfies Eq. (20) is a feasible 
solution of Eq. (19). 

A (D ehii K Ci + K T Ci D T ek ^ < -d < 0 (20) 

Vfc £ [h], Vi £ [n c ], 

where K Ci = diag{Kj,j £ c,;} and n c is the number of 
diagonal block for matrices in J. 




Fig. 2. IEEE 14 bus system. The dash lines are the com¬ 
munication links. Every generation bus, possibly with local 
load, is replaced by one VSI in this work for simplicity. 

Proof. From Eq. (20), there exist U = U T >- 0 such 
that 

-f- U Au i) < — d < 0, Vi G [Z]. 

By Schur compliment, Mi(e, £, £) R 0 if and only if 
{A\ hi U + UA lhl + £U) + e(UU + (I) R 0. (21) 


consider a general case where part of inverters may be 
disconnected from the microgrid abruptly due to some 
severe events. Let Vj f C V/ be a set of normal oper¬ 
ating inverter buses so that V/ \ Vi f is the set of dis¬ 
connected inverters. The voltage magnitude and phase 
angle dynamics at buses i g V/\ Vj f become unknown 
and are categorized as load buses. For this reason, we 
will partition the buses by V = VV ( U Vl,, for which 
Vi j contains all the normal operating inverters, while 
\>l s = V \ Vi f consists of all the other buses including 
load buses Vl or disconnected inverter buses V/ \ Vj f . 
Consider the case that the communication network re¬ 
mains intact when some inverters are disconnected from 
the microgrid. The communication graph between the 
operating inverters is defined as Qf := (Vi f ,£f), where 
£f = E c \{{i,j}\i £ V/\V/ / }. Let Lf be the Laplacian of 
Qf . The control law of the “fault” microgrid is reduced 
from Eq. (8) to 

Xi{t) = Ki Ejevi, L f(i,j)Si(x(t )) 

( 22 ) 

subj. to x(t) G X c , Vt > 0. 

Notice that due to the assumption on intact communi¬ 
cation network, the control law can be autonomously 
transformed to Eq. (22) in response to the change of mi¬ 
crogrid operating conditions. Let xi f := {xi,i € Vr f } 
and Si f := [S,,. i G Vi f }. The microgrid dynamics with 
the controller (22) is rewritten as follows 


We can find e, £ G R+ such that Eq. (21) holds when 
d G R + and the desired result follows. □ 


±i f (t) = K If L f Si f (x(t )) 
subj. to x(t) G X c , Vf > 0, 


(23) 


Lemma 4 allows us to find K by solving Eq. (20) in¬ 
stead of Eq. (19). The computational complexity is 
greatly reduced because the number of BMIs associated 
to Eq. (20) increases with respect to l = max i£ r„j U 
instead of n. In most microgrid networks, l is much 
smaller than n. Solving Eq. (20) instead of Eq. (19) 
enhances the practicability of the proposed distributed 
control method, especially for large-scale microgrids. 

4 Flexible Operation of Microgrids 

In this section, we study the plug-and-play feature of 
DERs under the proposed distributed control frame¬ 
work. The robustness under communication failure is 
also analyzed to solidify the proposed distributed con¬ 
troller. 

4-1 Plug-and-Play Feature 


where Kj f = diag{K.i,i G V/,}, and Lf = Lf (g) / 2 . 
Denoted Of as the proportional power sharing space 
of the inverters i G Vi r If Qf stays connected, the 
null space of Lf remains equivalent to Of. The “re¬ 
duced” microgrid (23) can therefore be analyzed through 
a similar way discussed in the last section. Define 'Df = 
CO{D eif , • ■ • , D eif ) as Jacobians of Si, when x G X c 

and l f as the number of vertices of Vf. The control ob¬ 
jectives (3)-(5) of the reduced microgrid follow if all the 
hypotheses in Corollary 1 hold except that Eq. (19) is 
replaced by the following 

maximize £> 

(24) 

subject to M i j{e f ,Cf,£f) R 0, 
where e/,£/ G M + and 
i 1MJ = [1,0]TjT 1 D gif Ki f LfTf[1, 0] T , 

M / x _ + u fAu,ij + tfCfl + tfUf U f 

lJ ’ L u f -tf 1 . 


The plug-and-play feature of the DERs refers to the 
property that one DER can be plugged or unplugged to a 
microgrid without re-engineering the entire control. We 
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Uf = Uj >- 0, i G [//]. 







The orthogonal matrix Tf of system (23) is defined sim¬ 
ilar to T as shown in Eq. (12), where the last two coordi¬ 
nates span the space of proportional power sharing. If K 
is found such that BMI constraints of Eq. (19) and (24) 
are both satisfied for V/ and all possible Vi f respectively, 
the proposed controller has the plug-and-play feature. 

In fact, we will show a lemma proving that the Eq. (24) 
is feasible if Eq. (19) is feasible and Qf is connected. 


The condition A (H r f) C R_ c is sufficient for the feasibil¬ 
ity of BMI constraints in Eq. (24) by Schur compliment 
and the property of Lf shown in Eq. (26). □ 


With Lemma 5, we can conclude stability of the “fault” 
microgrids by only solving optimization problem (19), 
which reduces the computational efforts to find a proper 
K for various microgrids operating conditions. 


Lemma 5 If Qf is connected and e*,K*,U*,( are 
feasible solution of Eq. (19), then by setting Ki f = 
diag{I\*,i G Vi f }, BMI constraints in Eq. (24) reduce 
to LMIs and are feasible. 


Proof. By definition of the Laplacian matrix, L and Lf 
are positive semidefinite. In addition, because Q c and Qf 
are connected, 


T~ 1 LT = 


Li 0 

0 0 


, 0 -< Li G K 2 («/- 1 )x2(m-i) ) (25) 


Tj'LfTf 


~Lfi 0 

0 0 




(26) 


4-2 Loss of Communication links 

We consider the case of communication failures in this 
subsection. Let Qi := (V/, £i) be a simple communica¬ 
tion graph when communication links £ c \ £i are failed. 
Denote L; as the Laplacian of Qi . We can write down the 
microgrid dynamics with the controller (8) by replacing 
L by Li. The analyses in the Section 3 apply and we can 
derive conditions to ensure objectives (3)-(5) under loss 
of communication links. 

Consider a more general scenario when both part of 
the DERs and communication links failed. Let Qfi = 
( Vi f ,£fi ) and Lfi be the graph and Laplacian under the 
operating condition, respectively. Assume the commu¬ 
nication links 


Since e*,K*, U* solve Eq. (19), the following condition 
holds for all i G [Z] 

[U*,0\T- 1 D g .K*T[L 1 ,0] T + (27) 

[L 1 ,C\T- 1 {K*) t D T e fT[U*,C\ T ~< -(eU*U* + CU*). 

The matrix inequalities are derived by applying Schur 
complement to BMI constraints of Eq. (19) and sym¬ 
metric property of L\. Due to Eq. (27), the eigenvalues 
of Hi := D gi K* + (K*) T Dl can be transformed into 
the following form by spectral decomposition 


where X(Hn) C R_ c , c > 0. The matrix Hu is a bijective 
mapping from H t to Hi, where Hi is some Euclidean 
space such that Hi 2 R 2n/ \ O. Notice that D eif Ki f + 

Kj f D ^ is a principal submatrix of Hi and 

H r ,i := [I,0}T^(D eiJ K If + Kf f Dlj)Tf[L0} r 

is a linear mapping from M 2 ! I \Of to R 2 I V/ jI \ Of, 
where M 2 I V/ ^ \ Of C Hi . As a result, H r , i is a princi¬ 
pal submatrix of Hu with a proper basis of H t . By the 
Cauchy interlace theorem, X(H r i ) C R_ c for all i G [If]. 


{(i, j),i G V If ,j€Vi\V If }C£\ £ fl (28) 

remains intact so that the inverters can still update the 
Laplacian Lfi autonomously. The desired properties in 
the scenario of loss of DERs and communication links 
tie closely to the feasibility of Eq. (24). We conclude this 
section by the following theorem 

Theorem 2 If the hypotheses in Corollary 1 hold, then 
there exists an non-empty X c , s C X c such that for some 
given loads and initial x(0) G X c s , the microgrid con¬ 
verges exponentially to the set where the control objec¬ 
tives including (3)-(5) are all satisfied if any one of the 
following operating conditions holds. 

(a) Microgrid operates normally with the dynamics (8), 

(b) The DERs i G V/\ Vi f are disconnected with the mi¬ 
crogrid, where Qf remains connected. The microgrid 
dynamics is described by (23). Voltage dynamics of 
buses V \ Vi f is bounded in the form of Proposition 1, 

(c) The communication network is reduced to a con¬ 
nected simple graph Qi. The microgrid dynamics is 
described by Eq. (8) with L replaced by L*i, 

(d) Both parts of the DERs are unplugged and communi¬ 
cation links are failed. The reduced communication 
graph Qfi stays connected. The dynamics is described 
by Eq. (23) withLf replaced by Lfi. Voltage dynam¬ 
ics of buses V \ Vi fl is bounded in the form of Propo¬ 
sition 1. 
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Proof. Every microgrid condition listed above has its 
own reduced order system in the form of Eq. (15). To 
prove the theorem, it is sufficient to show that all the 
reduced order systems are robustly stable with the hy¬ 
potheses of Corollary 1. Namely, the robust stability 
condition of every reduced order system in the form of 
Eq. (19) is feasible. The operating conditions (a) and (b) 
have their robust stability conditions hold due to Corol¬ 
lary 1 and Lemma 5. The operating condition (c) has 
the same dynamics to Eq. ( 8 ), only L is replaced by Z/. 
Notice that the eigenspace of 0 for both L and Li are the 
same. As a result, when L is replaced by L; in Eq. (19), 
the new BMI constraints remain feasible with the gain 
K optimizing Eq. (19). The robust stability conditions 
of operating scenario (c) follow. The feasibility of BMI 
constraints for operating scenario (d) follows if Eq. (24) 
is feasible because the only difference between the two 
is the replacement of Lf with Lfi, where similar argu¬ 
ments from scenario (a) to (c) apply. □ 

Remark 7 All the operating conditions listed in Theo¬ 
rem 2 require that the communication graph stays con¬ 
nected. Otherwise, the null space of the Laplacian matrix 
is bigger than the space of proportional power sharing. 
Proportional power sharing is only guaranteed between 
inverters which are connected by the remaining com¬ 
munication network. Toward this end, tertiary control 
may help establish proportional power sharing. The nor¬ 
malized power difference between “disconnected” com¬ 
ponents is first collected. The tertiary controller then 
broadcasts a new voltage reference to offset the power 
difference. Proportional power sharing can be gradually 
established with an aid of tertiary control. 

5 Simulation Results 

In this section, we validate the proposed controller by 
simulating IEEE 14-bus system shown in Fig. 2. The 
circuit specification can be found in [49]. Every gener¬ 
ator in IEEE 14-bus example is substituted by VSIs. 
We assume that the voltage magnitude deviation should 
be less than 6% of its nominal value for every bus. In 
addition, the branch angle difference limit is set to be 
less than 7 = 15 deg. We find a feasible solution K of 
Eq. (19) by solving Eq. (20) and maximizing 6 , where T> 
is replaced by J to reduce the complexity as discussed 
in Remark 6 . The calculated feedback gain K is shown 
in Table (1). Notice that a linear constraint on I\ is 
imposed 

— b m < KLs < 6m, Vs G t)(H), (29) 

where H represents active and reactive power capac¬ 
ity of inverter buses, 6m and b m are the upper and 
lower bounds of ramp of phase angle and voltage. We 
select 6 m and b m such that 49.7 < w < 50.3(Hz) and 
\E\ < 0.05(p.u/s) for all inverters. The constraint on fre¬ 
quency is standard for power networks [35]. The upper 


Table 1 

Feedback Gains 



Unit 

DER 1 

DER 2 

DER 3 

DER6 DER8 

Ki{ 1,1) 

mrad/s 

-5.6 

-1.8 

-10 

-8.8 -190.7 

Ki( 1,2) 

mrad/s 

10.2 

4.9 

8.6 

85 -24.6 

A, (2,1) 

mV/s 

-2.7 

-1 

-1.2 

-23.8 0 

Ki( 2,2) 

mV/s 

-0.1 

-1.4 

-10.8 

-16.2 -40 


bound on \E\ makes the voltage at inverter bus main¬ 
tain relatively static under any loading condition, while 
it gradually evolves to the steady state where the power 
sharing follows. Note that constraint (29) can be consid¬ 
ered as imposing an upper bound of | \K\ [ 2 - Adding con¬ 
straint (29) to Eq. (19) only changes the optimal value 
but not the feasibility. As a result, we can always select 
6 m and b m only based on physical requirements for the 
distributed controller design. 

We first simulate the case when all the inverters oper¬ 
ate normally with an abrupt change of load at bus 10 
at time t = 1. As shown in Fig. 3(a)- 3(e), the voltage 
magnitude always lies inside the desired range. In addi¬ 
tion, the mesh microgrid successfully reaches a satisfac¬ 
tory new steady state with proportional power sharing 
and synchronized frequencies. The second case that we 
consider is the failure of the inverter 1. From Theorem 2, 
K solves Eq. (19) can maintain the stability and achieve 
the requirements (3)-(5) when inverter 1 is disconnected. 
The simulation results are shown in Fig. 4(a)- 4(e). It 
can be seen that the rest of the inverters can carry over 
the original power injection from the inverter 1 while the 
desired properties are preserved. 


6 Conclusion 

In this paper, we have presented a distributed control 
method to coordinate VSIs in an island AC microgrid. 
Instead of using the conventional droop controller induc¬ 
ing steady state deviation of both frequency and volt¬ 
age magnitude, our design controls the microgrid to the 
steady state where frequency synchronization and pro¬ 
portional active and reactive power sharing hold while 
respecting the voltage regulation constraints. The suffi¬ 
cient conditions of the convergence to the steady state 
can be approximated into solving optimization problem 
with BMI constraints. The design is robust with respect 
to system parameter variations. Unlike most of the ex¬ 
isting droop based controller, our design can be applied 
to the lossy mesh microgrids. Relaxation of the commu¬ 
nication requirements and further reduction of the com¬ 
putational complexity are considered as important fu¬ 
ture works. Extension of the controller for higher order 
harmonic loads, unbalanced, and other types of nonlin¬ 
ear loads are also among our future works. 
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(a) Active power response of the inverters 
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(b) Reactive power response of the inverters 
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(b) Reactive power response of the inverters 
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(d) Voltage magnitude of buses 



(c) Frequency response of the inverters 
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(e) Phase angle difference between connected buses 

Fig. 3. Controlled microgrid response when the load at bus 
10 changes abruptly at time t = 1 


Time (s) 


(e) Phase angle difference between connected buses 

Fig. 4. Controlled microgrid response when the inverter 1 is 
failed and disconnected from the microgrid at 1 seconds. 


A Proof of Proposition 1 


Taking time derivative of Eq. (A.l) 


We first consider the case of constant impedance loads. 
Given the impedance at every load bus i, the power in¬ 
jection at bus i is determined by local voltage and is 
written as P*(xi ) and Q*{xi), for all i C Vl- Applying 
KCL at every load bus with power flow equation (1), we 
have 


0 = P*(x i ) + P i (x) 
0 = Q*(xi) + Qi(x) 


for all i € Vl- 


(A.l) 


®2til — fl,x.Xl T fL,xXLi (A.2) 

where fi tX and /z, iX denote the Jacobian matrices fol¬ 
low from Eq. (A.l). For any x, Hivlh < Ki^ll^/lb al¬ 
ways follows from Eq. (A.2) if /l, x is invertible. In the 
case that /l,x is not invertible, a change of coordinate 
is applied, let Xl = Tl, x xl, such that xl = [xj ,x^] T 

and Xl 2 span the zero eigenspace of /l,x- Let K 2 , x be 
the smallest positive eigenvalue of the Hermitian matrix 
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fI,xfL,x, then 

"• T 2, / • T rT £ • • X .r X ./* 

K2,xX Ll X Ll < X L f L>x f LtX X L = Xjf Ix f^ x X j 

= ^*| I't'Li 112 ^ ^3,x | |*t/1 12 j ^3,x G • 

Since there exist a simple linear mapping from to 
*l 2) I 112 < M.xIliiJb, we derive 

\\xl lb = 11*1,112 < | |*I/i 1 12 + ||iL 2 ||2 

< Kl,*||*j||2, « l,x G K+ 

for the case that //, )X is not invertible. For x € T c , we 
choose k = ma,x x( zx c «i,x and ||*l ||2 < /c||a:j||2 follows. 
The derivation also applies to constant power loads, only 
P* and Q* at the load buses become constants instead 
of the function of local voltage profile, which completes 
the proof. 
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